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Eigenfunctions of Id disordered Hamiltonian with constant 
imaginary vector potential are investigated. Even within the 
domain of complex eigenvalues the wave functions are shown 
to be strongly localized. However, this localization is of a very 
unusual kind. The logarithm of the wave function at differ- 
ent coordinates x fluctuates strongly (just like the position of 
Brownian particle fluctuates in time). After approaching its 
maximal value the logarithm decreases like the square root of 
the distance (In \ipmax/4>\) 2 ~ \x — xo\. The extension of the 
model to the quasi-ld case is also considered. 

PACS numbers: 72.15.Rn, 05.40.+j, 74.60.Ge 



The non-Hermitean disordered Hamiltonians with 
imaginary vector potential were introduced in Ref. |l|] 
for the description of vortex depinning from columnar 
defects by a transverse magnetic field in superconductors 
(later the analogous equations were used in the context 
of investigation of some problems of biological growth 
H). The features of this unusual Hamiltonians have at- 
tracted the immediate and wide interest P-^0[. From 
the theoretical point of view the most exciting was the 
prediction of the existence of a mobility edge in Id in 
such systems. However, the main attention in the ex- 
isting literature was paid to the solution of eigenvalue 
problem in the imaginary potential. The spectrum of 
this non-Hermitean Hamiltonian was shown to consist of 
two segments of real eigenvalues and two arcs of imag- 
inary eigenvalues (see insert on the Fig. 1). The do- 
main of real energy was calculated in the initial paper 
Ref. H . The analytic formula relating the real and imag- 
inary part of complex eigenvalues for Id model with weak 
disorder was found in Ref. || (some other analytic ap- 
proaches and approximations were considered in Refs. 

While all eigenstates with real energy turn out 
to be naturally localized, the complex eigenvalues are 
usually associated with the extended states. Neverthe- 
less, in this note we show that even after the transition to 
complex spectrum the wave functions remain sufficiently 
localized in Id imaginary vector potential. This localiza- 
tion, however, is of a very unusual kind. The logarithm 
of the wave function changes randomly with the increase 
of x. The maximum of the wave function corresponds 
to the maximal positive deviation of this random path 
from the initial position. After approaching maximum 
the logarithm of the wave function decreases also ran- 
domly being (ln|VWx| -ln|^(x)|) 2 ~ \x - x rnax \. We 
will call the localization of the eigenvectors of this type 



the Stochastic Localization. 

The Schrodinger equation for a particle in one channel 
disordered ring with constant imaginary vector potential 
has the form 
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V(x)ip = Ei(j 



(1) 



The boundary conditions are ip(x + L) = ip(x) (the 
closed ring). For wire with open ends the vector po- 
tential may be washed out by the gauge transformation 
if)(x) — > exp(hx)t/j(x). Our results do not depend on the 
concrete choice of the disordered potential V(x) (it may 
be white noise or something else). Also for analytic esti- 
mates we consider the case of weak disorder. This means 
that both the wave length and h~ x are small compared 
to the mean free path. 

We will consider only the states with large complex 
energy. Eigenvectors with small real energy may be ob- 
tained by gauge transformation from the usual localized 
states with h = 0. The most clear difference of the Hamil- 
tonian Eq. (|J) from the usual Id disordered case is that at 
h large enough the left going and right going waves are 
no more divergent. The corresponding energies for the 
momentum ±fc are E± = k 2 — h 2 ± 2ikh. As a result one 
may look for the solution of Eq. ([!]) of the form (analog 
of the WKB method in Hermitean Quantum Mechanics 

H) 



ip — exp(jcr) 



(2) 



The small corrections due to back scattering may be es- 
timated by the method which we use below for the inves- 
tigation of multichannel case. The derivative expansion 
for a now gives 



a = cr + <7i 



(3) 



cr = / WE- V - ih)dx' ; cri = - \n{E - V) 



The weak disorder may mix only the plain-wave states 
with momentum k' close to k. Therefore, only slow gar- 
monics of the potential Vkk' came into play, which jus- 
tifies the use of WKB approximation . For our pur- 
poses it will be enough to consider only the first term 
(To of the expansion Eq. (||). The quantization condi- 
tion now takes the form a(L) — <t(0) = 2im, which in the 
first order in V leads to the energy E = (k + ih) 2 + (V) 
with k = 2-Kn/L and (V) = L^ 1 J Q Vdx. The function 
u(x) is almost real. However, due to a weak disorder it 
acquires a small fluctuating imaginary part. Due to this 
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imaginary part the logarithm of the modulus of the wave 
function (again in the leading order in V) behaves like 



In 



h 2 



{V{x') ~ (V))dx' 



(4) 



Thus we see that our wave function may be large (expo- 
nentially!) or small depending on the occasional sign of 
the integral. We see also that the momentum fc emerges 
in the Eq. (||) only in the overall factor (fc 2 + h 2 )^ 1 . As 
a result all wave functions turn out to be localized at 
the same place. This feature of localization in the strong 
imaginary potential is in sharp contrast with what hap- 
pens in the Hermitean Id disordered Hamiltonians, where 
eigenfunctions with close energy strongly repel each other 
in the coordinate space. In the case of white noise disor- 
der V(x)V{x') — DS(x — x') one may find how fast the 
logarithm Eq. (0) changes from point to point. Namely 



(ln|^)|-ln|^(t/)|) 5 



h 2 D\x-y\(L 



y\) 



A(k 2 + h 2 ) 2 L 



(5) 



In particular this equation allows us to estimate, how the 
wave function decreases after approaching the maximum. 
One may introduce via the Eq. (H) the typical size of the 
wave packet (analog of the localization length) 



2\2 
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(fc 2 + h 2 ) 

h 2 D 



(6) 



This localization length formally is of the same order of 
magnitude as the usual Anderson localization length in 
Id disordered wire £4 ~ 1/D. However, the nature of 
this localization is completely different. In Anderson 
case all wave functions decrease at large distances like 
exp(— x/£a). In our case the wave functions decrease like 
exp(— const yx/^h) and even the constant in the argu- 
ment of the exponent fluctuates with x and from sample 
to sample (no thermodynamic limit). 

Our results Eqs. (|||^) are valid only far from the tran- 
sition from real to complex spectrum in the limit of weak 
disorder. In order to study the more general situation, 
numerical simulations were performed. For numerical 
computations it is natural to consider the tight binding 
variant of the model ([l]) 



t (e h ip n+ i + eS/>n-i) + V„ip n = Eip„ 



(J) 



The approximate solution of this equation may be easily 
found for the case /i> 1 (compare also with |l(J) 



E k « -t\e h - i<t ' + e- h+i4 



-te h -^ ; = 2tt- (8) 



= exp I in<i> + e-' l +^ ^ (V m - (V)) + 0{e~ 2h ) \ . 

I m=0 ) 

Here integer N and fc are the size of the closed 
chain and the number of solution (momentum), (V) = 



N^ 1 Sn=o Vn- F° r concrete calculation we have used 
N = 300, t = 1, h = 0.4, and the random potential uni- 
formly distributed within the range —1.5 < V n < 1.5. 
Thus we consider the case of intermediate or even strong 
disorder. The eigenvalues on the complex plane for 
one realization of disorder are plotted on the insert on 
Fig. 1. The left segment of real eigenvalues includes 29 
low energy states. The more interesting for us are the 
squared wave functions |V , n' ) | 2 shown on the same figure 
(the wave functions are normalized to unity in the usual 



1 and only the values \ip n \ > 10 5 are 



shown). We have shown the few typical states well below 
the phase transition i = 9, 10, 11 and the states around 
and above the transition i = 27, 30, 36, 50. The solutions 
are numbered in accordance with increase of ReE. First, 
one sees that all the states below and above the "mobil- 
ity edge" are clearly localized. Second, this localization 
has a very different form below and above the transition. 
The consecutive real eigenfunctions i — 9, 10, 11 practi- 
cally do not overlap. Alternatively, the states above the 
transition (and even slightly below) turns out to be very 
similar. This features of the wave functions are in explicit 
accordance with what we expected from the analytic es- 
timates for weak disorder. On the Fig. 2 we have shown 
the In \ip { n ] | 2 for i = 4, 50, 150. As we have mentioned be- 
fore, up to exponentially small corrections the real wave 
functions may be obtained from the solutions with h = 
via ip n (h) = e nh ifj n (h = 0). In accordance with this pre- 
diction we see on the Fig. 3 for i = 4 that the logarithm of 
the wave function decreases linearly from the maximum 
with different slopes to the left and to the right. Also 
as we expected, the logarithm of the complex functions 
resembles to large extent the Brownian path. Moreover, 
due to the Eq. (j|) at the middle of the zone {4> — tt/4) 
the contribution of the leading order in V changes only 
a phase of the wave function but not an amplitude. The 
Eq. (@) is valid only for h ^S> 1 and this peculiarity at 
the center of the band should not necessary survive in 
our simulations. However, as we see, the effect of lo- 
calization is strongly suppressed (but is not washed out 
completely!) for i = 150. 

In order to introduce the analog of orthonormal basis 
in non-Hcrmitean case one has to consider two sets of 
eigenvectors: Hij)R = Ei/jr and ipjJ'H = Eip^. Up 
to now in fact we have considered only the ipR. For 
our specific choice of the Hamiltonian Eq. ([!]) evidently 
ipL(E,h) = iPr(E* ,—h). Moreover, in the leading ap- 
proximation, with only cto left in the Eq. (^|), one simply 
has i/jl{x)iPr(x) = const. Thus we see that the growth 
(exponential!) ofipR at the maximum is compensated by 
the decrease of ipL. This compensation is not an artifact 
of approximation Eqs. (|^||). For example, following the 
authors of the Ref. |Q, one may introduce the conserv- 
ing current J = rp' L ipR — iPl^'r + IhipLipR dJ/dx = 0. 
This means that the strong effect of stochastic localiza- 
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tion may be washed out if the observables are bilinear 
in ipn and ipL, but should be seen if the observables de- 
pend on only one tpR or ipi,. Surprisingly one finds the 
examples of both kinds within the most popular appli- 
cation of non-Hermitean Quantum Mechanics. Namely, 
the authors of Rcfs. |l],[l5|] have used the Hamiltonian 
Eq. (|l|) for description of Abrikosov vortex depinning by 
transverse magnetic field in the hollow superconducting 
cylinder with random columnar defects. The probabil- 
ity to find vortex at the point x on the transverse slice 
r of the cylinder with length L T was shown to have the 
form P{x,t) = Z- 1 {^f \e-^- T ^ H \x){x\e~ TH . For 
very large L T only the ground state contribution survives 
in the evolution operator e~ LrH . In this case naturally 
0,1! P(x,t) ~ ^ R {x)^ L {x) for r,L T -r ~ L T . We see 
that the effect of localization above depinning transition, 
which we consider in this paper, could hardly be seen 
inside the cylinder fl7|| . However, at the top of cylinder 
P(x,L T ) ~ ipn(x) and at the bottom P(x, 0) ~ i}>l{x) 
the "stochastic localization", which we have introduced, 
turns out to be of 100% importance. We leave the more 
detailed discussion of the localization above depinning in 
case of many vortices for a separate publication Jl8| . 

Up to now we have studied the effect of imaginary vec- 
tor potential only in the Id systems. Consider now the 
closed strip of the width I and length L with random 
weak disorder V(x, y)V(x', y') — D 2 5{x — x')S(y — y') 
and imaginary potential along the ring (the Hamiltoni- 
ans of the kind Eq. (p]) in more than Id have been also 
considered in Refs. The Hamiltonian for the 

strip differs from that of the Eq. ([!]) by the only addi- 
tional term —d 2 /dy 2 . It is convenient to look for the 
solutions of this Schrodinger equation of the form 

i>(x,y) = ^2ip n (x)yf2/isui(q n y) , (9) 

with q n — mc/l. The analog of the Eq. (Rl) now reads 



(d/dx- 



(10) 



where the longitudinal energy s n =_E-q 2 n . In the case of 
white noise disorder V nn (x)V nn {x') — (3D 2 /2l)8(x — x') 
and for to and n different V nm (x)V mn (x') = (D 2 /l)8(x — 
x'). At least for small I (or large e n — e m ) and for disorder 
weak enough one may consider the hopping V nm in Eq. 
( |l0| ) as a perturbation. The zeroth approximation for 
the wave function in this simple case is again — 
exp{ J (V e « — Km — ih)dx'} with some given n and E = 
£ n+Qn = (k+ih) 2 +q 2 . Physically interesting, however, is 
the range of validity of this single chanel approximation. 
The corrections to ipn^ are described by the equation 



[- {d/dx - h)* 



(1) _ 



-V ^ (0) 



(ii) 



The solution of this equation may be found by the 
Green's function method. Let for definiteness e m — e„ > 



(note that Ime m = Ime n ). In this case the Green's func- 
tion for the Eq. ( |TT| ) G(x — y) equals zero at x > y and 
has a cusp at x = y. The solution reads 



oo X 1 (x-y) 



X 2 — Ai 



V mn (y)^\y)dy 



(12) 



where Ai 1=1 2h — ik and A2 ~ ik + (e m — e n )/{h — ik) 
(k is related with E as we have shown above and we 
suppose that \e m — e n \ <C £„). For |n — m\ <C n one has 
Re\ 2 <C k,h. However, if (k 2 + h 2 )\e m — e„\ > hD 2 /l 
(compare with Eq. (|^)) the amplitude of the exponent 
|exp(A2£)| still vary with x much faster than ^^(ie)!- 
In this case the averaged value 



Do 



1 



32nhq n - 



WW 



(13) 



The analogous formula for e m — e„ < is obtained after 
the replacement n — m — ► \n — m\. In the similar way 
one may find the correction to the direct potential V nn 
induced by the hopping V nm in the second order of per- 
turbation theory. This correction in its turn renormalizes 
the speed of stochastic growth of the amplitude of wave 
function Eq. (||) and the localization length Eq. (||). The 
generalization of the Eq. (j^) for the strip gives 



(ln|V>i%)M%)l) 2 = 



(14) 



y\D 2 



Do 



4(fe 2 + h 2 )l 



E 



32nhq '—f \n 



Here \x — y\ <SC L. The sum over m in the r.h.s. is log- 
arithmically divergent. The summation should be cut at 
I n — to I ~ kl (we suppose that k ~ q ~ h). Thus we see 
that at small I the growth of the wave function simply 
reproduces that of the single chanel case Eq. (j^) with the 
effective localization length proportional to the width of 
the strip £/, ~ /. Only at D 2 \n(kl)/hq ~ 1 the features 
of wave functions in a strip became sufficiently different 
from those in 1-chanel wire. One may interpret this re- 
sults as an indication of the existence of exponentially 
large localization length on 2d plane with imaginary vec- 
tor potential. 

One may easily repeat the calculations Eqs. (f^lif) also 
in the case of thick 3d wire. The result of such general- 
ization will be the formula analogous to the Eq. (|l4|) and 
£,h ~ l 2 /D^ (I is the typical transverse size of the wire). 
Also like in 2d case the new physics starts only at ex- 
ponentially thick wire D 3 hi(kl)/q ~ 1. However, the 
concrete behaviour of the perturbed wave function may 
be completely different for 2c? and 3d. In a strip due to a 
complex energy E = (kii + ih) 2 + q'j_ there is no degener- 
acy of different transverse channels. Thus the only new 
phenomena we may expect with increase of the width of 
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the strip is localization of wave functions in transverse 
direction due to the mixing of channels with very close 
q±. In 3c? case the analog of the Eq. ([l3]) describes the 
diffusive mixing of many channels with close (q±) but 
completely different q±. 

In summary, we have shown in this paper that the tran- 
sition from real to complex spectrum in Id disordered 
systems with imaginary vector potential is not a dereal- 
ization transition. However, the nature of localization be- 
low and above the transition is completely different. The 
most interesting is the localization at complex energies, 
there the wave function decreases (increases) like the ex- 
ponent of y 'x/^h and this exponent even has no a well de- 
fined thermodynamic limit. As for physical applications, 
the localization at complex energies becomes invisible for 
observable values bilinear in left- and right- eigenvectors 
of non-Hcrmitcan Hamiltonian, but is of 100% impor- 
tance for observables depending on only one of them ip^ 
or ipR. For example, the effect which we consider will 
lead to the strong modification of the distribution of flux 
lines at the ends of superconducting cylindrical shell with 
random columnar defects. For the 2d(3d) generalizations 
of the model the same effect takes place with the effec- 
tive localization length £h proportional to the transverse 
size(area) at least until the system became exponentially 
thick. 
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FIG. 1. The IV'n'l 2 for the model Eq. (7) for states with 
numbers i = 9, 10, 11, 27, 30, 36, 50. We use N = 300, t = 1, 
h = 0.4, and —1.5 < V n < 1.5. The first vector with complex 
energy is i — 30. All states are clearly localized. However, 
the states far below the "mobility edge" differ strongly from 
each other, while the complex states look quite similar. Insert 
shows the eigenvalues on the complex plane. 
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FIG. 2. The In \ip^ \ 2 for i = 4 ("localized") and i = 50, 150 
("stochastic") states. The data for i = 150 was multiplied by 
the factor 5. The Brownian curves for complex states and 
different left and right slope for i — 4 are clearly seen. 
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